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ABSTRACT 

We report the presence of high significance diffuse radio emission from the Triangulum Aus¬ 
tralis cluster using observations made with the KAT-7 telescope and propose that this emis¬ 
sion is a giant radio halo. We compare the radio power from this proposed halo with X-ray 
and SZ measurements and demonstrate that it is consistent with the established scaling rela¬ 
tions for cluster haloes. By combining the X-ray and SZ data we calculate the ratio of non- 
thermal to thermal electron pressure within Triangulum Australis to be X = 0.658 ±0.054. 
We use this ratio to constrain the maximum magnetic field strength within the halo region to 
be 5max,haio = 33.08/iG and compare this with the minimum field strength from equiparti- 
tion of 5min,haio = 0.77(1 to place limits on the range of allowed magnetic field 

strength within this cluster. We compare these values to those for more well-studied systems 
and discuss these results in the context of equipartition of non-thermal energy densities within 
clusters of galaxies. 

Key words: galaxies: clusters: intracluster medium - galaxies: clusters: individual: Triangu¬ 
lum Australis 


1 INTRODUCTION 

A number of galaxy clusters are sources of diffuse radio emission 
that can be classified as either radio halos or radio relics (e.g., Fer- 
etti et al. 2012, and references therein). The radio emission is syn¬ 
chrotron radiation produced by relativistic electrons with Lorentz 
factors of the order of 10^ that move in juG magnetic fields. 

Giant radio halos have sizes of 1 — 2 Mpc, are located at the 
centres of clusters, have fairly steep spectra and are not usually 
observed to have significant polarization (e.g., Feretti et al. 2001; 
Bacchi et al. 2003). Synchrotron emission from such large volumes 
requires that local particle acceleration is effective throughout the 
cluster (Jaffe 1977). Although the basic observational properties 
of radio halos have been established (e.g., Feretti et al. 2012), the 
formation mechanism of radio halos is still unclear (e.g., Brunetti 
et al. 2008; Donnert et al. 2010a,b; Macario et al. 2010; Brown 
& Rudnick 2011; EnBlin et al. 2011; Brunetti et al. 2012; Zan- 
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danel et al. 2012; Arlen et al. 2012). Theories that explain their ori¬ 
gins include primary models, in which an existing electron popula¬ 
tion is re-accelerated by turbulence caused by recent cluster merg¬ 
ers (Brunetti et al. 2001; Petrosian 2001), and secondary models, 
in which relativistic electrons are continuously injected into the 
ICM by inelastic collisions between cosmic rays and thermal ions 
(e.g., Dennison 1980; Blasi & Colafrancesco 1999; Dolag & EnBlin 
2000; Miniati et al. 2001; Keshet & Loeb 2010). Combinations of 
both acceleration mechanisms have also been considered (Brunetti 
& Blasi 2005; Dolag et al. 2008; Brunetti & Lazarian 2011). 

Since few radio telescopes cover the very low declinations, 
most radio halos are found in the Northern sky. The only radio 
halo known below a declination of —40 deg is the bullet cluster 
(Liang et al. 2000.) In order to extend the sample of radio halos, 
we started from the BAX [^cluster catalog, selecting those objects 
with declination < — 40 deg, T > 4 keV, z < 0.5, and some evidence 
of a merger either from the ROSAT images or from the literature. 

^ http://bax.ast.obs-mip.fr 
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Figure 1. 1328 MHz KAT-7 image of the Triangulum Australis region. Left: KAT-7 data are shown as greyscale and contours for the full field without point 
source subtraction. The half power point of the KAT-7 primary beam is shown as a circle and the synthesized beam is shown in the bottom left corner with 
dimensions of 3.67 x 3.41 arcmin. Greyscale data are saturated at 30mJybeam“^ in order to highlight the low surface brightness diffuse emission. Contours 
are shown in increments of 5 CTnns from 5 Right: KAT-7 data are shown as contours, with intervals as in the left-hand figure; SUMSS data are shown as 
greyscale, saturated at 400mJybeam“^ in order to highlight low surface brightness emission. The half power point of the KAT-7 primary beam is shown as a 
circle. In both maps the KAT-7 G^ms = 1.84mJybeam“^ and no correction is applied for the KAT-7 primary beam response. 
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Figure 2. 1328 MHz point source subtracted KAT-7 image of the proposed 
halo emission within Triangulum Australis. Contours are shown in incre¬ 
ments of 1 CTrms from 4 CTj-ins, where Grms = 1.84mJybeam“^ No correction 
for the primary beam response has been applied to these data. 


This resulted in a sample of eight clusters. In this paper we present 
the first of these: Triangulum Australis. The Triangulum Australis 
cluster is a relatively nearby (z = 0.051) bright, hot system, which 
was overlooked in the optical band due to its low Galactic latitude. 
It was first discovered as an X-ray source (McHardy et al. 1981). 
The Triangulum Australis cluster has been observed with XMM- 
Newton (60 ks, Markevitch et al. 1996) and it was found that this 
cluster has a hot (12 keV) core at its centre that is most likely pro¬ 
duced by a merger. 

Finally, the cluster is close enough (z = 0.051, 5' ~ 300 kpc) 


that even at low resolution a radio halo could be resolved. As part of 
the development of MeerKAT (Booth et al. 2009), a scientific test 
array, the Karoo Array Telescope (KAT-7), has been constructed 
and commissioned at the same site. In this paper we report the dis¬ 
covery of a giant radio halo with the KAT-7 array, showing the po¬ 
tential of the array to image extended, low-surface brightness ob¬ 
jects. 

Throughout this paper we assume a ACDM cosmology with 
Ho = 67.3 kms“^ Mpc“^ Llm = 0.32, and = 0.68. All images 
are in the J2000 coordinate system and all errors are quoted at 1 a. 


2 OBSERVATIONS 

The KAT-7 telescope consists of seven 12 m diameter dishes, 
equipped with cryogenically cooled receivers working between 
1.3 GHz and 1.8 GHz with an observational bandwidth of 
256 MHz. The dish distribution is optimized for a Gaussian UV dis¬ 
tribution, with highest weighting given to the optimisation param¬ 
eters of 4-hour tracks at 60 degrees declination (de Villiers 2007). 
The maximum baseline separation is 192m and minimum spacing 
is 24m. 

Triangulum Australis was observed as part of general com¬ 
missioning for the KAT-7 instrument four times between 2013 Feb 
and 2013 June at a central frequency of 1.328 GHz, giving a total 
integration time of approximately 40 hours. For each observation, 
primary calibration was performed using PKS 1934-638, while sec¬ 
ondary gain calibration used PKS 1718-649. 


2.1 Data reduction 

The native KAT-7 data comes in the Hierarchical Data Format 
(hdf5). Once converted into measurement set (ms) format using 
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in-house software, the data were reduced using the CASA pack- 
ag^ Channels contaminated by known RFI were flagged immedi¬ 
ately, followed by automated flagging using the CASA RFLAG rou¬ 
tine, looking at both auto- and cross-polarisation components. After 
flagging, the data were calibrated following standard practice. Flux 
densities were set using PKS 1934-638, tied to the Perley-Butler- 
2010 flux density scale in SETJY. 

MS-MFS deconvolution was carried out using the CLEAN task 
in CASA over a 2° x 2° fleld-of-view (FOV; 1.5 times the FWHM 
of the KAT-7 primary beam). Imaging was performed by initially 
using a mask based on sources from the SUMSS catalogue (Mauch 
et al. 2003) with 843 MHz flux densities exceeding 15 mJy, before 
removing the mask to allow deconvolution of the whole held. The 
resulting Stokes I image is shown in Fig. [^and has an rms noise of 

= 1.84mJybeam“^, which is measured using the rms in the 
central region of the source subtracted image. These data are confu¬ 
sion limited at the resolution of KAT-7. Predictions of the expected 
confusion level for KAT-7 at this frequency are slightly lower than 
the measured rms noise in these data, cJconf — 1.4mJybeam“^ 
(Riseley et al. 2014). We attribute this difference to the enhanced 
source population towards galaxy clusters, relative to the held. 


3 RESULTS 

Diffuse emission towards the Triangulum Australis cluster is visi¬ 
ble in the KAT-7 data at a signiflcance of > lOcr over an extent of 
several arcminutes and a major axis of approximately 1 Mpc within 
the 5 (7 contour, see Fig. The diffuse radio emission is coinci¬ 
dent with the X-ray emission towards this cluster, although a slight 
offset 2 arcmin) exists between the peaks of the radio and X- 
ray emission. The KAT-7 image has astrometry for the radio point 
source population in this held consistent with previous high resolu¬ 
tion surveys (SUMSS; Mauch et al. 2003), see Fig.[^(right), which 
suggests that this offset is due to the differing nature of the physi¬ 
cal processes responsible for the X-ray and radio emission: X-ray 
emission predominantly traces the density of the thermal gas pop¬ 
ulation within the cluster, whereas radio emission traces the cos¬ 
mic ray electron population and magnetic held strength distribu¬ 
tion. Offsets between the peak surface brightness of different emis¬ 
sion mechanisms in disturbed clusters are not uncommon, see e.g. 
Rodriguez-Gonzalvez et al. (2011). We do not make a further phys¬ 
ical interpretation of this offset in this work. 

From X-ray studies, it has previously been proposed that Tri¬ 
angulum Australis is a merging system due to its high central tem¬ 
perature (Markevitch et al. 1996) and therefore likely to host a giant 
radio halo. We propose here that the diffuse radio emission detected 
with KAT-7 is associated with that halo. It is possible that the exten¬ 
sion of the radio emission seen towards the North of the cluster may 
be due in part to an unresolved cluster relic; however, given the low 
signiflcance of this protrusion we do not try to separate these fea¬ 
tures. We note that neither diffuse nor compact emission is present 
in SUMSS (Mauch et al. 2003) data towards the proposed halo 
emission. This reduces the possibility that the emission detected 
with KAT-7 is due to a collection of unresolved point sources. 


WWW.nrao.edu/casapy 


3.1 Point source removal 

A large number of point sources are also detected within the KAT- 
7 held of view, see Fig. In order to reduce any confusing ef¬ 
fect on the diffuse emission identifled with the halo, these sources 
were used to solve for direction dependent calibration solutions be¬ 
fore being subtracted from the visibility dataset. At the frequencies 
of our observations, the cross-correlation between voltages from 
pairs of antenna are affected by a series of moderate but complex 
baseline-time-frequency Direction Dependent Effects (DDE). They 
might include atmospheric effects, pointing errors or dish deforma¬ 
tion. 

A large variety of solvers have been developed to tackle these 
kinds of calibration issues. Here we do not attempt to physically 
characterize the DDEs but instead use a Jones-based solver. This 
type of solver constitutes the most widely used family of algo¬ 
rithms for direction dependent calibration, and aims at estimating 
the apparent net product of the various effects mentioned above. 
Recently, algorithms have been developed (see e.g. Yatawatta et al. 
2008; Noordam et al. 2010), that estimate a Jones matrix per time- 
frequency bin per antenna, per direction. The well known prob¬ 
lems associated with this type of technique are (i) ill-conditioning 
and (ii) computational cost, both being due to the larger number 
of degrees of freedom used to solve the problem (compared to the 
direction-independent case). The first of these issues can affect the 
scientific signal by suppressing un-modelled flux, while the cubic 
computational cost with the number of degrees of freedom can put 
strong limitations on the affordable number of direction-dependent 
parameters. 

The Jones-based solver utilized here (Tasse et al. 2014) is a 
DDE variant of the StefCal approach (Salvini et al. 2014). It op¬ 
erates using the concept of iteratively solving for linear systems in 
a similar manner to traditional non-linear least-squares solvers, but 
by using an alternative iteration scheme, significantly improving 
convergence speed and robustness. 

Eor the data presented here, in order to increase the signal in 
each direction, we clustered the sources in 5 direction-based groups 
by using a Voronoi tessellation algorithm and computed a scalar 
direction dependent Jones matrix every 15 minutes. We verified that 
this strategy was not driving suppression of the unmodelled flux by 
using incomplete sky models. 

Eollowing direction-dependent calibration, point sources 
above a signiflcance of 7.5 cJrms were then subtracted directly from 
the visibility data. Point source subtracted data were then re-imaged 
using natural weighting in order to enhance the signal-to-noise of 
the low surface brightness halo on large scales. The point-source- 
subtracted image is shown in Eig.[^ 

Sources with flux densities above 7.5 Oi-ms within the 
50 per cent power point of the KAT-7 primary beam are listed in 
Table where they are cross-referenced with the v2.1 SUMSS 
catalogue (Mauch et al. 2008; hereafter M08). Errors on KAT- 
7 flux densities as liste d in Table [T] are calculated as cr = 
\l + ^fit Where multiple SUMSS sources are as¬ 

sociated with a single KAT-7 detection, due to the large difference 
in resolution between these two surveys, the combined flux density 
is used to calculate the spectral index value; in this case, uncertain¬ 
ties on the SUMSS data points are combined in quadrature. The 
average spectral index for the sources in this list is ^ 343 ^^^ = 
0.66 ± 0.43, typical of optically thin non-thermal emission. 
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Table 1. Sources detected at a significance of ^ 7.5 CTrms within the half power point of the KAT-7 primary beam with SUMSS counterparts. Column [1] lists 
a numerical designation for each source; column [2] lists the fitted Right Ascension, with the fit error on this position listed in column [3]; column [4] lists 
the fitted Declination, with the fit error on this position listed in column [5]; column [6] lists the fitted peak flux density for each source from the KAT-7 data; 
where sources are unresolved this value is listed as column [7] lists the integrated flux density for each source from the KAT-7 data; column [8] lists the 
SUMSS flux density for coincident sources; column [9] lists the SUMSS source designation for the listed SUMSS flux densities where names are composed of 
the truncated SUMSS J2000 coordinates; column [10] lists the measured spectral index between the SUMSS and KAT-7 flux densities; column [11] identifies 
notes on the KAT-7 data fitting. 


No. 

RA 


Dec 


‘5peak,1328 

‘5int,1328 

‘^int,843 

SUMSS 

^1328 

^843 

Notes 


[J2000] 

[s] 

[J2000] 

[arcsec] 

[mJy bm“^] 

[mJy] 

[mJy] 

Identifier 



001 

16 33 13.2 

±0.31 

-64 23 00 

±4.69 

47.4±3.8 

82.9±6.2 j 

f 34.7±1.5 
[ 50.9±3.1 

J163309-642322 ] 
J163317-642157 J 

^ 0.07±0.19 

(i;2) 

002 

16 33 25.8 

±0.42 

-64 28 22 

±4.33 

77.5±5.2 

116.5±7.6 

131.6±4.1 

J163327-642832 

0.27±0.16 

(2) 

003 

16 33 28.7 

±0.54 

-64 07 46 

±8.19 

27.8±2.1 

38.4±2.9 

60.8±2.1 

J163331-640805 

1.01±0.18 

(2) 

004 

16 34 37.7 

±0.38 

-64 41 12 

±5.72 

27.6±2.5 

54.9±3.7 

47.8±1.8 

J163434-644040 

-0.30±0.17 

(2) 

005 

16 34 53.4 

±0.52 

-64 14 25 

±7.79 

- 

19.1±2.5 

40.8±1.6 

J163453-641420 

1.67±0.30 


006 

16 35 01.1 

±0.42 

-63 58 44 

±6.32 

98.5±5.8 

117.7±6.9 j 

f 152.4±6.0 
[ 29.9 ±2.9 

J163457-635838 j 
J163523-640040 J 

^ 0.96±0.15 

(2) 

007 

16 36 30.0 

±0.20 

-64 35 22 

±2.96 

_ 

54.1±3.6 

90.8±2.9 

1163629-643515 

1.14±0.16 


008 

16 36 51.7 

±0.25 

-65 08 07 

±3.75 

66.1±4.4 

82.4±5.4 

101.3±3.2 

1163652-650808 

0.45±0.16 


009 

16 38 10.5 

±0.44 

-65 04 29 

±6.64 

198.0±13.3 

290.0±19.6 j 

f 347.9±13.9 
[ 47.0±1.9 

1163808-650409 ] 
1163751-650414 J 

^ 0.68±0.18 

(3) 

010 

16 38 13.1 

±0.36 

-63 55 19 

±5.44 

16.6±1.0 

27.5±2.3 

32.6±1.4 

1163816-635536 

0.37 ±0.21 

(3) 

on 

16 38 31.8 

±0.19 

-64 41 04 

±2.90 

89.0±5.1 

92.6±5.3 

105.0±5.5 

1163830-644043 

0.28±0.17 


012 

16 39 07.5 

±0.25 

-65 07 20 

±3.79 

203.8±12.6 

337.0±20.8 j 

f 339.6±13.6 
[ 95.1±6.1 

1163913-650804 ] 
1163903-650513 J 

^ 0.56±0.16 

(3) 

013 

16 39 24.2 

±0.10 

-64 05 13 

±1.50 

64.7±3.8 

68.1±4.0 

89.7±4.5 

1163924-640520 

0.61 ±0.17 


014 

16 40 05.5 

±0.23 

-64 26 42 

±3.43 

79.2±5.9 

124.5±8.9 j 

f 47.6±3.0 
[ 96.4±5.7 

1164000-642639 ] 
1164007-642717 J 

^ 0.32±0.19 

(3;4) 

015 

16 40 24.3 

±0.96 

-65 05 35 

±14.43 

_ 

32.2±3.0 

55.9±2.1 

1164033-650529 

1.21 ±0.22 


016 

16 41 45.9 

±0.05 

-64 34 02 

±0.77 

296.2±15.1 

302.4±15.4 

443.5±13.4 

1164145-643407 

0.84±0.13 


017 

16 42 38.4 

±0.36 

-64 20 30 

±5.38 

- 

26.9±2.6 

53.6±1.9 

1164239-642050 

1.51 ±0.23 


018 

16 43 55.8 

±0.11 

-64 40 16 

±1.63 

148.8±8.2 

154.8±8.5 

238.9±7.3 

1164354-644019 

0.95±0.14 



Note 1: Adjacent source SUMSS J163254-643254 may also contribute emission at 843 MHz. Fitting may be affected. 

Note 2: Diffuse component evident in KAT-7 data. 

Note 3: Closely adjacent source. Fitting may be affected. 

Note 4: Additional component listed in original SUMSS catalogue (Mauch et al. 2003) but not in later version (Mauch et al. 2008). 


3.2 Flux density estimation 

Integrated flux density measurements for the proposed halo were 
made using the source subtracted images in order to avoid con¬ 
tamination from the point source population. The region of dif¬ 
fuse emission that we associate here with the radio halo of Tri¬ 
angulum Australis is extended and irregular. The centroid of the 
diffuse emission is located at J 16h38m48.5s —64°20'13" and the 
peak of the diffuse radio emission at J 16h38m52.5s —64°22'0l", 
see Figure Within the 3 oi-ms contour the diffuse emission has 
dimensions of 850 x 990 arcsec (East-West by North-South) and a 
major axis of 1022 arcsec, where 1 arcsec is 1.035 kpc at z = 0.051. 
The CASA task IMFIT applied to this target returns a value for 
the integrated flux density of Si^fn = 186 ± 15 mJy. However, this 
method involves fitting a Gaussian to our target region, which is 
significantly non-Gaussian in morphology. Consequently we also 
use an aperture photometry technique to extract the integrated flux 


density using the FITFLUX code (Green 2007). This method fits 
a tilted plane to the edges of a defined aperture before subtract¬ 
ing this plane and integrating the remaining flux density. Using 
this method on data corrected for a Gaussian primary beam with 
a FWHM of 1.31 degrees, we And an integrated flux density for the 
radio halo of ^gtflux = 130±4mJy, where the error on the fitted 
value is calculated using the standard deviation of the recovered 
flux density from multiple apertures of varying dimension. We cal¬ 
culate our complete uncertainty on the integrated flux density as 

a = ^+ (0.055i„t,fitfiux)^ to give a final integrated 
flux density of 5'int,halo = 130 ± 8 mJy. 


4 SCALING RELATIONS 

The power-law relationship between radio power, Pj-ad^ ^nd X-ray 
luminosity. Lx, is well known for clusters hosting haloes (e.g. Gio- 
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vannini et al. 2000; Feretti et al. 2012) and is commonly charac¬ 
terised as Pj-ad where d has values of approximately 1.5 — 2.1 
(Brunetti et al. 2009). A further power-law relationship is also 
known linking radio power and the integrated Compton-y param¬ 
eter, Fsz, determined from observations of the Sunyaev-Zel’dovich 
(SZ) effect. Unlike the X-ray luminosity, which depends on the 
propertes of the thermal components within the cluster, Fsz is pro¬ 
portional to the total electron pressure integrated along the line of 
sight (Colafrancesco et al. 2003). Consequently, the correlation of 
radio power to integrated Compton-y is of particular interest as it in¬ 
dicates the relationship between the non-thermal electron pressure 
component (characterised by Pj-ad _^^l5(a+l)/2v-(«+')/2 - 

^non-thermal^ electron pressure. Following 

Colafrancesco et al. (2014), we denote the ratio of these quantities 
as X, where 

Aotal ~ thermal 3“^non—thermal ~ 3“-^)thermal* (1) 

As can be seen from Figure the radio power from Trian¬ 
gulum Australis is consistent with the known scaling relations. 
Here we use the sample of Colafrancesco et al. (2014), which 
extends the sample of Basu et al. (2009). Integrated Compton-y 
values are taken from the Planck catalogue (Planck Collaboration 
2013), which measures the cylindrical volume integrated Compton- 
y, Fcyi = FszP^A’ within an aperture of P = 5P5oo- At this radius the 
cylindrical integrated quantity is equivalent to the spherically inte¬ 
grated quantity, Fsph (Arnaud et al. 2010). Furthermore, it can then 
be trivially related to Fr 5 oo as Fr 5 oo = F 5 R 500 x ^(l)/^(5), where 
7(1) = 0.6552 and 7(5) = 1.1885 (Appendix 2 of Arnaud et al. 
2010 ). 


4.1 Non-thermal pressure fraction 


The value of the ratio X can be determined from the X-ray luminos¬ 
ity, Lx and the integrated Compton-y value, as calculated at Psoo^ 
the radius at which the cluster density profile is equal to 500 times 
the critical density of the Universe, Pcrit(^) = 377^ (z)/8 ttG, where 
77 (z) = 77oL (z) and L (z) = (1 + z) ^ ^ . Again following 

Colafrancesco et al. (2014), this relationship is. 


^sph,R500^(^)^^"^ = 


,5/4 


r5/4 


where the constants Fq and Lq may be found from 


( 2 ) 


Fo = ^^^^Cpmp500pcritWeO,gV'i(A) (3) 

/ 271 \ 9 9 

Lo = 4;rC2 ( —GMmp500perit j «?o.g^ (4) 


where p = 1.14 is the mean molecular weight, G is the gravita¬ 
tional constant, mp is the proton mass, oy is the Thomson scattering 
cross-section and C 2 has the value 1.728 x lO'^^Ws”^ 

(Rybicki & Lightman 1985) with 



This calculation assumes that the global cluster density profile is 
modelled by a /3 -model, with index /3 and Tq = AR 500 and therefore 
has a central electron number density of 


_ 3/3/B500pcrit 

2 A 2 pemp 


( 7 ) 




Figure 3. Cluster halo scaling relations. Top: X-ray luminosity and radio 
power scaling relation; Bottom: Integrated Compton-y and radio power 
scaling relation. Data are taken from Colafrancesco et al. (2014; black 
points) with the exception of Triangulum Australis (red squares), which 
has properties as determined in this work. The fitted power-law relations 
are taken from Colafrancesco et al. (2014) and shown as dashed lines. 


where /b is the baryon fraction, here assumed to have the value 
/b = 0.175 (Planck Collaboration 2013). 

For Triangulum Australis we use /3 = 0.63 ± 0.02, rc = 3.5 ± 
0.2arcmin (Markevitch et al. 1996) and A = 0.3, consistent with 
Colafrancesco et al. (2014). Combining these with the X-ray lumi¬ 
nosity and integrated Compton-y, we find that X = 0.658 ± 0.054. 


4.2 Maximum magnetic field strength 

The thermal electron pressure of the cluster within R 500 is ex¬ 
pressed as 

Ah,500 = 500^3^500 (8) 

= ^+X)-‘. (9) 

Consequently, the average non-thermal pressure within R 500 can be 
calculated as Lnon-th,500 = ^Ah,500- This additional pressure con¬ 
tribution to the SZ effect will come from the non-thermal particle 
population, with other kinetic forms of non-thermal pressure such 
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as tubulence and bulk motions contributing to the kinetic SZ (kSZ) 
effect. Turbulence is generally assumed to be the dominant form of 
non-thermal pressure (Vazza et al. 2012), but due to the directional 
nature of the kSZ effect and the multiple line of sight reversals ex¬ 
pected for a turbulent medium, the net turbulent contribution to the 
kSZ effect is expected to be negligible, as is that of bulk motions 
(Sunyaev et al. 2003). 

Furthermore, it is expected that magnetic pressure is sub¬ 
dominant and will not be greater than non-thermal particle pres¬ 
sure (e.g. Lagana et al. 2010; Brunetti & Jones 2014). Under this 
assumption, one may calculate an upper limit on the strength of the 
cluster magnetic field such that 

^max,500 ^ x/SttP^oii— th,5005 (10) 

where B is the magnetic field strength in Gauss and pressure is 
measured in Barye. From Eq. using a representative temper¬ 
ature of 7500 = lOkeV (Markevitch et al. 1996) and noting that 
^e,500 = 500/bPc(z) /this provides an upper limit on the av¬ 
erage magnetic field strength of (P)max,500 = 14.50pG. Alterna¬ 
tively, using the integrated Compton-F value and Eq. this pro¬ 
vides an upper limit on the average magnetic field strength within 
P 500 for Triangulum Australis of (P)max,500 = 19.46pG. 

Since the magnetic field strength is expected to vary as a func¬ 
tion of cluster radius we convert our value of (P) max,500 to be more 
representative of the field within the halo region, which has a radius 
of 5 arcmin, less than 50% Rsm- In order to do this we assume that 
the dependence of magnetic field strength on density follows 

(11) 

V ^^6,0 / 

(Bonafede et al. 2013) where rj =0.5, consistent with that deter¬ 
mined for the Coma cluster. Combining this dependence with the 
standard /3-model formalism for the radial density distribution, this 
gives 

B(r)oc (^l + (r/rcfj . (12) 

The magnetic field strength determined from Eq. is a volume- 
averaged quantity within Psoo- Assuming spherical symmetry 
and (P)max,500 = 19.46pG, the maximum average magnetic field 
strength within the halo radius, r^, is given by (P)max,halo = 
33.08 jUG. 


4.3 Minimum magnetic field strength 


The minimum magnetic field strength, P^iin in Tesla, within the 
halo region can be calculated, assuming equipartition, from its ra¬ 
dio power following 


3/Uo G(a)(l+k)Prad 
2 Vf 


2/7 


(13) 


(Longair 2011) where k is the ratio between the energy of heavy 
particles (protons) and the electrons, / is the filling factor used to 
describe the fraction of the volume, V, occupied by radio emitting 
material and G(a) is defined as 


G{a) = 


a{p)(p-2) 

(7.4126 X 10-‘®)-(7’-2) 


r ^,-( p - 2)/2 _ -{ p - 2 )/ 2 l { p -\)/2 

L^min ^max J v' 


2.344 X 10 


1-25 


(1.253 xl0^^)“(^“i)/^ 


( 14 ) 


with 

'/^r(p/4+19/12)r(p/4-l/12)r(p/4 + 5/4) 

2 (p+l)r(p/4 + 7/4) ’ ^ ^ 

where p = I — 2a. Since we only have a measurement of the ra¬ 
dio power at a single frequency, we must assume a spectral in¬ 
dex, a. Here we use a = —1.5, consistent with previously mea¬ 
sured halo indices (Eeretti et al. 2012). This gives p = 4.0 and 
a{p) = 0.034; using Vmin = 10 MHz and Vmax = 100 GHz we find 
G{a) = 0.256v^-^. Eor our measured fiux density of 5i. 33 GHz = 
130 ± 8mJy, given the redshift of z = 0.051, the radio power is 
Prad = 0-6 X 10^^ WHz“^. Eor clusters of galaxies a value of /: = 0 
or ^ = 1 is typically used (Beck & Krause 2005); however, the same 
authors also propose that in fact a larger value of k{k» 1; such 
that n^/riQ ^ 100 — 300) is preferred by current models of cosmic 
ray production in galaxy clusters. Recent constraints using a com¬ 
bination of radio data and upper limits from gamma-ray observa¬ 
tions has shown that in galaxy clusters, n^/riQ is likely to be signif¬ 
icantly less than 100 (Vazza & Briiggen 2014; see also Guo, Sironi 
& Narayan 2014). Here we model the radio halo as a solid sphere 
with radius, = 5 arcmin. We assume that the volume is filled uni¬ 
formly and completely, such that / = 1. We find that, given these 
assumptions, P rnin = 0.77(1 +/:)^/^jUG. 

We note that although equipartition and minimum energy ar¬ 
guments are frequently used, they are subject to a number of strong 
assumptions. One of the strongest assumptions is the value of the 
parameter k, the effect of which we have explicitly factored in our 
estimates. A further issue is the strong local dependence of the ra¬ 
dio emissivity on magnetic field strength, which can cause Pmin to 
over-estimate the volume-average field strength for inhomogeneous 
magnetic fields. A more complete discussion of these assumptions 
is presented in Beck & Krause (2005). 


5 DISCUSSION 

Erom considering the combination of X-ray, SZ and radio data we 
are able to place both lower and upper limits on the magnetic field 
strength in the halo region. The upper limit in this case assumes 
that the magnetic pressure will not be greater than all non-thermal 
particle pressures. Our derived upper limit for the magnetic field in 
the halo region of Triangulum Australis is (5)niax,halo = 33.08 pG. 
Since the cooling time (via synchrotron emission) of electrons is a 
function of magnetic field strength, 

where 5 cmb — 3(1 + z)^ pG is the energy density equivalent mag¬ 
netic field strength of the CMB, this field strength would imply a 
short synchrotron cooling time of ~7 Myr. If the relativistic elec¬ 
trons of the radio halo are re-accelerated by first-order Eermi pro¬ 
cesses (as in Brunetti et al. 2001), the short loss time implies that 
the electron distribution evolves quickly, leading to an expedient 
decrease of the break energy of the electron spectrum. The short 
synchrotron loss time indicates that whatever mechanism powers 
the radio halo must still be active. Eor hadronic models, a short 
cooling time means that the radio emission must follow the ICM 
density distribution quite closely. This could potentially be probed 
by radio observations at higher angular resolution. Eor models of 
turbulent re-acceleration, it implies that the turbulence should have 
a high filling factor and be efficient in electron re-acceleration. 
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Here, low-frequency observations would be useful (see e.g. Figs. 7 
and 9 in Brunetti et al. 2001). 

For comparison we consider the Coma cluster of galaxies 
(A1656), which is a particularly well studied system with ex¬ 
tensive ancillary data available. Using /3-model parameters from 
Briel, Henry & Bohringer (1992), such that /3 = 0.75 ±0.03 and 
Tc = 10.5 ± 0.6 arcmin, an X-ray luminosity of (10.44 ±0.28) x 
10^^ergs“^ (Reichert et al. 2011), an integrated Compton-y of 
(0.1173 ± 0.0054) arcmin^ (Planck Collaboration 2011), a 1.4 GHz 
radio halo power of (0.72 ±0.06) x 10^"^ WHz“^ (Brunetti 2009), a 
halo radius of = 21 arcmin (Venturi et al. 1990) and assumptions 
consistent with those outlined above we find 5niin,eq = 0.46(1 ± 

llG, (5)max,R500 ~ 10.81/iG and (5)max,halo ~ 16.32/iG. Wc 
note that the value of X = 0.322 for Coma (Colafrancesco et al. 
2014) is approximately twice the vaue of dp/p = 0.15 value deter¬ 
mined by Fusco-Femiano et al. (2013) for Coma. These values are 
consistent given that p^/p ^0.5. 

Our equipartition value for the Coma cluster is consistent with 
that of Thierbach et al. (2003), who find 5niin,eq = 0-68 pG (with 
k = 1), allowing for varying cosmologies. Although it has been ar¬ 
gued that this value is an underestimate, due to the choice of k\ Beck 
& Krause (2005) suggest that this field strength could be as high as 
4pG, assuming rip/nQ = 1000. Magnetic field strength measure¬ 
ments for Coma have also been made using the Faraday rotation 
of polarized emission from its galactic population (Bonafede et al. 
2013). Unlike minimum energy equipartition measurements, Fara¬ 
day rotation, (j), provides a direct measure of the magnetic field 
strength along the line of sight (l.o.s.) such that (j) oc 
For the Coma cluster, Faraday rotation measurements indicate that 
the average l.o.s. magnetic field strength within Tc is 5|| = 4.7 pG 
(Bonafede et al. 2010), consistent with the limits set here using 
equipartition and the non-thermal pressure fraction. 

We note that the calculations outlined here and in Section |4] 
assume that the clusters under examination are well described by a 
/3-model. This assumption creates limitations in the situation where 
either the cluster gas density or non-thermal halo gas population de¬ 
viates significantly from spherical symmetry. In this situation the 
estimates for the non-thermal pressure (and hence the magnetic 
field strength) are likely to differ. There are two potential causes for 
such a situation in this case: firstly, this is a merging system and the 
assumed beta model is may not be a good representation; secondly, 
a possible unresolved radio relic could bias the radio power high. 
Given the low radio luminosity of the proposed Triangulum Aus¬ 
tralis halo relative to the general scaling relations, see Fig. the 
second of these scenarios seems unlikely; however, observations 
at higher resolution with improved sensitivity relative to currently 
available data are required to examine this possibility in more de¬ 
tail. The former scenario is likely to affect the results presented 
here, but to what extent is currently unclear. Further development 
of the methodology used to calculate the non-thermal fraction will 
be necessary to assess the impact. 


6 CONCLUSIONS 

We have used new observations with the KAT-7 telescope to make 
the first detection of a diffuse radio halo in the Triangulum Aus¬ 
tralis cluster. By combining these new radio data with complemen¬ 
tary data in the X-ray and SZ regimes, we have demonstrated that 
this cluster is consistent with the established scaling relations for 
clusters hosting haloes. In addition we have: 


• Used a combination of X-ray and SZ data to determine the 
ratio of non-thermal to thermal pressure within the cluster, which 
we determine to be X = 0.658 ± 0.054. 

• From this ratio of pressures we were able to determine 
an upper limit on the average magnetic field strength within 
^ 500 ’ (^)max,500 = 19.46jUG, and hence within the halo region, 
(^)max,halo ~ 33.08/iG. 

• We have compared these values with the lower limit equipar¬ 
tition value determined from the radio power, under stated as¬ 
sumptions, which we determine to be = 0.77(1 ±k)^/^;UG. 
Hence providing both lower and upper limits on the possible field 
strengths within the cluster halo region. 

• We use the well-studied Coma radio halo to contextualize 
these results and demonstrate that the range of values we calcu¬ 
late for the allowable magnetic field strengths are consistent with 
measurements made using alternative methods. 
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